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Sensitivity studies to the CP-violating parameters of the decay B® — > <j><f> with the LHCb exper- 
iment are presented. The decay proceeds via a b — > sss gluonic-penguin quark transition, which 
is sensitive to contributions from beyond the Standard Model particles. A time-dependent angular 
analysis of simulated data leads to an expected statistical uncertainty of 6° on any new physics 
induced CP-violating phase for a sample corresponding to 2 fb _1 of integrated luminosity. The 
expected precision on sin 2(3 from the related decay B° — > 4>Kg is also discussed. 
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I. INTRODUCTION 

The amplitude for the decay B® — > (jxj) is dominated 
by gluonic-penguin quark transitions b — > sss (Fig. [!}. 
Gluonic-penguin processes are sensitive to beyond the 
Standard Model particles that contribute within the loop. 
The e + e~ P-factories have measured sin 2/3 in nine P° 
gluonic penguin decay modes, such as B° — > 4>Kg and 
B° —> v'Kg [lj. All the measurements of sin 2(3 from 
these modes have values below that measured in b — > 
ccs transitions, but no individual measurement shows a 
significant deviation. 

The decay B° s — > 0</> is predicted to have a CP- 
violating phase less than 1° within the SM [2]. The 
dependence on Vt s in both the mixing and decay am- 
plitudes leads to a cancellation of the P°-mixing phase. 
Therefore, if any significant CP-violation is measured 
in P° — > 4>4> decays it is an unambiguous signature of 
new physics. The decay is of a pseudoscalar meson to 
two vector mesons, which leads to the final state being 
a CP-even and CP-odd admixture. Therefore, a time- 
dependent angular analysis is required to extract the CP- 
violating parameters of the decay. 

The paper is organized as follows. SectionHTlcontains a 
brief description of the LHCb experiment. The predicted 
event yields and background estimations are described 




FIG. 1: The main diagram contributing to the decay 



in Section Hill The CP sensitivity study is presented in 
Section Hvl The LHCb prospects with the related mode 
P° — » <j>Kg are discussed in Section [V] The conclusions 
are given in Section [VII 



II. THE LHCb EXPERIMENT 

The Large Hadron Collider (LHC) collides protons 
at a centre-of-mass energy of 14 TeV. The LHC pro- 
duces 10 12 bb quark pairs per nominal year of data-taking 
(10 7 s) when operating at an instantaneous luminosity of 
2 x 10 32 cm^s" 1 . 1 The LHCb spectrometer (3,0] instru- 
ments one forward region about the pp collision point. 
The forward geometry captures approximately one-third 
of all B hadrons produced and increases the probability 
of both B hadrons from the bb pairs being within the 
acceptance, which improves the efficiency of flavour tag- 
ging. 

A silicon vertex detector, with sensors perpendicular to 
the beam axis, is situated close to the interaction region 
in a secondary vacuum. The detector provides accurate 
determination of primary and secondary vertices leading 
to a proper-time resolutions of approximately 40 fs in 
hadronic P-decays such as P° — > <jxj). Tracking stations 
either side of a 1.2 T dipole magnet produce momentum 
measurements with an accuracy of a few parts per mille. 
There are two Ring-Imaging Cerenkov detectors, with 3 
different radiators, that allow identification of K ± from 
7T^~ over the momentum range 1 to 100 GeV/c. 

In addition, the detector includes an electromagnetic 
calorimeter, a hadron calorimeter and a muon detector. 
These components are critical for identifying large trans- 
verse momentum, px, electrons, photons, hadrons and 
muons from P-hadron decay in the initial hardware stage 
(Level-0) of the LHCb trigger. The Level-0 trigger re- 
duces the 40 MHz collision rate to 1 MHz. All data is 
then transferred from the detector to a dedicated CPU 
farm where the Higher Level Trigger (HLT) algorithms 
are performed. Initially an association between the high 
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This luminosity optimises the number of single interactions per 
bunch crossing. 
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FIG. 2: The B® mass of candidates reconstructed in the signal 
(blue) and inclusive bb simulation samples (red) before trigger 
selections. The normalisations are arbitrary. 



Pt objects that satisfy the Level-0 trigger and tracks with 
large impact parameter to the primary vertex is sought. 
If this is successful, exclusive and inclusive HLT algo- 
rithms are executed resulting in a 2 kHz output rate to 
disk. 



III. EVENT SELECTION AND BACKGROUND 
ESTIMATION 

The estimation of B® — > (jxf> yields and background 
has been performed on simulated data. The simulation 
of the pp collisions and subsequent hadronisation is per- 
formed by the PYTHIA generator The decay of any 
B hadrons produced is simulated by EVTGEN [(J and the 
detector response is performed by GEANT4 [7j|. The re- 
sulting data are processed by the complete LHCb recon- 
struction software. A dedicated sample containing events 
with a B® — ► <fxj) decay is used to evaluate the selection 
efficiency. Minimum bias events are selected rarely by 
the Level-0 trigger and the HLT. Therefore, an inclusive 
sample of 34 million events containing bb quark pairs is 
used to estimate the background. Due to the very large 
number of bb pairs that will be produced at LHCb, the 
inclusive bb sample only corresponds to approximately 
15 minutes of data taking at the nominal instantaneous 
luminosity; this leads to large uncertainties on the back- 



ground estimates. 

The selection process begins with the identification of 
4> candidates reconstructed from two oppositely charged 
kaons. Particle identification and px requirements are 
placed on the K and they must be consistent with pro- 
duction at a common vertex. The mass of the K + K~ 
must be within 20 MeV/c 2 of the <f> mass; the mass in- 
terval corresponds to approximately ±3cr m ^, where u m<j) 
is the mass resolution. 

B-meson candidates are reconstructed in events with 
two or more (j> candidates with a pr > 1.2 GeV/c 2 . The 
P-meson candidates must be consistent with production 
at a common vertex and this vertex must be well sepa- 
rated from the primary vertex. The B® candidate mass 
distribution is shown in Fig. [5] for signal and background 
samples. Events within 40 MeV/c 2 of the B® mass are 
considered as signal; the mass interval corresponds to 
approximately ±3er mB , where <J mB is the mass resolu- 
tion. Once Level-0 and HLT trigger selections have been 
applied there are 4000 signal events expected in every 
2 fb _1 of integrated luminosity, which corresponds to 
one nominal year of LHC operation. The event yield 
is calculated assuming the measured branching fraction 
B(B° S -> (jxji) = (U±l(stat.) ± 6(syst.)) x 10" 6 @]; the 
measured value lies within theoretically predicted range 
[i| H(ij |. The background remaining in the bb inclusive 
simulation sample is found to consist of combinatoric B® 
candidates. The background-to-signal ratio is bounded 
to lie between 0.4 to 2.1 at the 90% confidence level, 
with a central value of 0.9. 



IV. B" 



CP SENSITIVITY 



The magnitude of any new physics induced CP-phase, 
4>np, is extracted from a time-dependent analysis of 
the differential cross section with respect to the three 
transversity angles defined in Fig [3] The amplitude for 
the decay can be written in terms of three helicity am- 
plitudes H\(t) where A = 0, ±1. The helicity amplitudes 
are related to the transversity basis by A (t) = H (t), 
A||(t) = ^(P +1 + and A ± (t) = ±(H +1 - H^). 

The amplitudes Ao and An are CP even and A± is CP 
odd. The differential cross section is then given by: 



, ~, AV f\ T - | J 4„(i)| 2 /i(x^i^2)|A||(i)| 2 /2(x^i^2) + | J 4 ± (t)| 2 / 3 (x^i^2) 

d\d cos u\d cos U2 

+ S(A* (t)A ± (t))f 4 ( X , 01,0a) + R(^(t)A||(«))/ 6 (x, O u 2 ) + 9(>4fj (t)A x (t))fe(x, 0i,9 2 ) (1), 

I 

where fi (i = 1 — 6) are even angular functions as re- tors of these angular functions are sensitive to any CP- 
quired by Bose symmetry (Tlj . The time-dependent fac- violating phase. In principle the CP- violating phase can 



3 




FIG. 3: A schematic of the definition of the transversity angles 
61, 02 and x- 
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FIG. 4: The variation of the acceptance as a function of the 
proper-time, r. The acceptance function fit to the simulated 
data is e(r) = 

be different between the three transversity amplitudes. 
However, to simplify the analysis it has been assumed 
to be equal. Significant fine-tuning of the phases would 
be required for no effects to be observed if new physics 
induced phases are present. The time-dependent terms 
also depend on the strong phase differences 5\\q between 
A(t)± and A(t)\\ t o, the relative magnitudes of the three 
amplitudes and the mass (lifetime) differences between 
the B® mass eigenstates, Ara s (AF S ). 

Simulated data are generated to follow the differential 
distribution given in Eqn. [T]with the value of <Pnp cho- 
sen to be 0.2 rad. The strong phases are assumed to be 
<5|| = and 5o = tt; these values are motivated by naive 
factorization [12] . The magnitudes of the transversity 
amplitudes are set to the values measured in the anal- 
ogous channel for B° decays B° -> K*6 [H, Q3]. The 
value of Am s is taken to be 17 ps -1 [15[ and Ar s /r is 
set to be 0.15, compatible with current experimental con- 
straints [l6| and theoretical expectations \Yj\ . The signal 
sample size has a mean of 4000 events corresponding to 
an integrated luminosity of 2 fb _1 . 

The following experimental effects arc also simulated. 

• The background is assumed to be at level 90% 
of the signal with flat transversity angle and mass 
distributions, and an exponential lifetime distribu- 
tions. 



• The tagging power, e(l — 2w), where e is the tag- 
ging efficiency and to is the mistag rate, is assumed 
to be 9% which has been found in simulation stud- 
ies of other J5° hadronic decays [l8| . This is signif- 
icantly better than that for B° modes because the 
kaon associated with the B® hadronisation is also 
used. 

• The proper-time acceptance measured from the 
signal simulation sample is shown in Fig. |4j The 
reduced acceptance for short lifetimes is the result 
of trigger and selection requirements on the impact 
parameters of the B daughters. 

• The proper time and B° mass resolutions are 

estimated to be 40 fs and 12 MeV/c 2 , respectively. 
These resolutions are estimated from the signal 
simulation sample used for the selection studies. 

• The angular acceptance and resolution are as- 
sumed to be flat and to have negligible effect, re- 
spectively; this assumption is motivated by the 
studies of related channel B° s — > Jtptp [Ti| . 



Prob 


1 


Constant 


20.59+ 1.373 


Mean 


0.1941 ±0.006382 


Sigma 


0.1002^0.004403 




FIG. 5: The distribution of fitted value of (/>np for 500 simu- 
lated B° — > (fxf) experiments. 

Five-hundred samples of signal and background events 
are generated and a maximum likelihood fit is performed 
on each one. The parameters 4>nPi $\\, So and the mag- 
nitude of the transversity amplitudes are extracted from 
the fit; all other parameters are fixed. The distribution 
of the fitted value of 4>np for these 500 experiments is 
shown in Fig The average error on 4>np is 0.1 rad 
(5.7°). The pull distribution of the true value subtracted 
from the fitted value divided by the uncertainty is nor- 
mal. 

Sets of 500 experiments are produced varying the input 
parameters assumed. The variation of the uncertainty 
on cj)NP as a function of the B(B° — > (fxf), signal-to- 
background ratio and Ar s /T s is given in Table |TJ The 
variation of the assumed B leads to the expected statisti- 
cal scaling of the uncertainty. The uncertainty on 4>np is 
not degraded significantly until the background-to-signal 
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TABLE I: The variation of the uncertainty on 4>np as a 
function of the branching fraction, background-to-signal ra- 
tio (b/s) and Aryrv 
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ratio is greater than three. Increasing the value of Ar s /r 
leads to a reduction in the uncertainty because enhanced 
interference from the lifetime difference among the am- 
plitudes increases the sensitivity to 4>np- The values of 
4>np, the proper-time resolution and the relative magni- 
tudes of the transversity amplitudes are also varied; all 
these had negligible effect on the sensitivity to 4>np- 

V. B° -> <j>K% CP SENSITIVITY 

Simulation studies of the decay B° — » 4>K S have also 
been performed. The expected yield per 2 fb _ of in- 
tegrated luminosity is 800 events. These yields do not 
include Kg daughters without measurements in the sili- 



con vertex detector, approximately two-thirds of K s from 
these decays, which are not reconstructed in the current 
HLT algorithms. Algorithms to perform this reconstruc- 
tion are currently being developed. The background-to- 
signal ratio is estimated to be 2.4 from the bb inclusive 
simulation. 

A time-dependent analysis of the B° — > 4>K S is re- 
quired to extract the sensitivity to sin 2/3. A toy sim- 
ulation study is performed to extract the sensitivity to 
sin 2/3. The tagging power is assumed to be 5% and the 
proper-time resolution is taken to be 60 fs from the simu- 
lated signal sample. A 10% K + K~ 5*-wave contribution 
is also included in the fit. The uncertainty on sin 2(3 is 
expected to be 0.32 for a data sample corresponding to 
2 fb _1 of integrated luminosity. 



VI. CONCLUSIONS 

Sensitivity studies to CP- violation in the decay — > 
(jxp have been presented. A sample of data corresponding 
to an integrated luminosity of 2 fb _1 gives an uncertainty 
of 6° on any new physics induced CP phase. Varying the 
assumptions used within reasonable ranges changes the 
predicted statistical uncertainty between 4° and 13°. The 
largest statistical uncertainty results from decreasing the 
branching fraction by a factor of four. The measurement 
of sin 2(3 from the B° — > 4>K S has also been investigated. 
The sensitivity is expected to be 0.32 with a data set 
corresponding to 2 fb _1 of integrated luminosity; this is 
of the same order as the current sensitivity of the e + e~ 
B-factories [2fj|. |2lj. Therefore, in conclusion, B® — > 4>4> 
is the most sensitive mode with which to study gluonic- 
penguin B decays with LHCb. 
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